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We  obtain  a  query  lower  bound  for  quantum  algorithms  solving  the  phase  estimation  problem. 
Our  analysis  generalizes  existing  lower  bound  approaches  to  the  case  where  the  oracle  Q  is  given  by 
controlled  powers  Qp  of  Q ,  as  it  is  for  example  in  Shor’s  order  finding  algorithm.  In  this  setting  we 
will  prove  a  Q(log  1/e)  lower  bound  for  the  number  of  applications  of  QP1  ,QP2, ....  This  bound  is 
tight  due  to  a  matching  upper  bound.  We  obtain  the  lower  bound  using  a  new  technique  based  on 
frequency  analysis. 
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I.  INTRODUCTION 


We  study  lower  bounds  for  the  phase  estimation  prob¬ 
lem.  In  this  problem  we  are  given  a  unitary  transfor¬ 
mation  Q  as  a  black-box  and  we  know  that  |g)  is  an 
eigenvector  of  Q,  i.e. 

Q\q)=e2^\q),  p  €  [0,  1).  (1) 


We  want  to  determine  the  phase  p  up  to  precision  e. 

The  quantum  phase  estimation  algorithm  approxi¬ 
mates  p  given  1^)  and  is  the  main  building  block  in 
Shor’s  factoring  algorithm,  the  counting  algorithm,  and 
the  eigenvalue  estimation  algorithm  [1  7]. 

The  main  element  of  this  algorithm  are  controlled  pow¬ 
ers  of  Q,  which  we  define  as  follows.  Let  Q  be  a  t  qubit 
unitary  transformation  and  |  tp)  an  arbitrary  t  qubit  state. 
For  l  =  1, . . . ,  c,  and  p  £  N  we  define  the  c+t  qubit  trans¬ 
formation 


Wf{Q)  \x\...xc)  \ip) 


\x\...xc)  \ip)  Xi  =  0 

\xi...xc)  Q'p  \ip)  xi  =  1 


We  call  Wf{Q)  a  (controlled)  power  query.  If  the  trans¬ 
formation  Q  is  clear  from  the  context,  we  will  just  write 
Wf  =  Wf(Q). 

This  notation  allows  us  to  write  the  phase  estimation 
algorithm  in  a  compact  form,  see  figure  1.  The  algorithm 
returns  an  approximation  p  of  the  phase  p  of  |  q). 
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FIG.  1:  The  quantum  phase  estimation  algorithm  in  power 
query  notation.  H  is  the  Hadamard  gate,  Wp  a  power  query 
as  in  equation  (2),  and  T~t  is  the  inverse  quantum  Fourier 
transform  on  T  qubits. 
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It  is  well  known  that  T  =  loge-1)  power  queries 
suffice  to  approximate  p  up  to  e.  In  this  paper  we  study 
whether  it  is  possible  to  improve  on  the  performance  of 
the  phase  estimation  procedure,  i.e.,  we  ask  what  is  the 
minimal  number  of  applications  of  WF  to  estimate  p  up 
to  e. 

Theorem  1.  Any  quantum  algorithm  estimating  the 
phase  p  of  an  eigenvector  |g)  of  matrices  Q  up  to  pre¬ 
cision  e,  with  Q  from  the  class 

Q\q),t  =  {Q  ■  Q  is  a  unitary  t  qubit  transform, 

|g)  is  an  eigenvector  of  Q}.  (3) 

has  to  use  fl(logi)  power  queries. 

We  prove  Theorem  1  in  section  IV. 

The  query  cost  of  the  phase  estimation  algorithm  may 
be  given  by  counting  each  application  of  W^Q)  as  p 
applications  of  Q ,  i.e.  the  query  cost  of  this  algorithm  is 

1  +  2  +  4  +  .  •  •  +  2t_1  =  2t  —  1. 

For  certain  problems,  like  order-finding,  it  is  possible 
to  exploit  some  knowledge  about  Q.  Here  Q  \y)  = 

| xy  mod  N)  for  a  certain  fixed  x,  and  therefore 

Q2  \y)  =  | x2J V  m°d  A")  =  |  (x2J  ^  y  mod  IV) 

is  easy  to  compute  by  repeated  squaring  and  modular 
multiplication,  see  e.g.  [2],  In  this  case  we  can  use  power 
queries  IT/2  with  essentially  the  same  cost  as  an  applica¬ 
tion  of  Q.  Thus  we  can  execute  the  phase  estimation  al¬ 
gorithm  with  query  cost  T  and  have  exponential  speedup 
for  the  query  cost:  from  2T  —  1  to  T.  Let  us  stress  that 
this  speedup  only  applies  if  the  cost  for  computing  Q 23 
is  similar  to  that  for  computing  Q. 


II.  PRIOR  WORK 

Quantum  query  complexity  has  been  important  in 
quantum  computing  since  Grover’s  search  algorithm, 
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which  is  provably  superior  to  classical  algorithms  [8,  9] 
in  the  number  of  queries.  The  first  lower  bound  result 
was  given  in  [10],  which  used  an  adversary  argument. 

Our  lower  bound  approach  is  based  on  the  ideas  of  the 
“polynomial  approach”  of  Beals  et.  al.,  [11-14],  Other 
approaches  include  the  quantum  adversary  argument  and 
its  generalizations  [15-17]. 

These  approaches  only  cover  problems  concerning  Boo¬ 
lean  functions,  so  we  have  to  extend  them  to  numerical 
problems.  This  has  been  done  through  extensions  of  the 
polynomial  method  [18,  19]  and  has  been  applied  widely 
to  integration  [20-22],  path  integration  [23],  approxima¬ 
tion  [24-26],  and  eigenvalue  estimation  [6,  7]. 

In  this  paper  we  apply  the  approach  of  [19]  to  the  phase 
estimation  problem.  Instead  of  using  a  maximum  degree 
argument,  which  is  not  applicable  in  the  case  of  arbitrary 
powers,  we  will  develop  a  new  lower  bound  technique 
based  on  frequency  analysis. 


IV.  GENERAL  CONTROLLED  ARBITRARY 
POWER  QUERIES 

We  consider  arbitrary  powers  p±,  . . .,  pr-  This  requires 
us  to  introduce  a  new  proof  technique.  To  illustrate  the 
idea  consider  the  phase  estimation  algorithm,  performed 
as  in  figure  1  with  c  =  T  =  3  control  qubits. 

0  I)W?W2W^(H®3  0  I)\0)  | q) . 

Let  us  trace  through  each  of  the  steps  in  this  algorithm 
(we  neglect  normalization  factors). 

1.  (|0>  +  |1)  +  |2)  +  |3)  +  . . .  +  |7))  | q) 

2.  (|0)  +  e2”^  |1)  +  |2)  +  e2'Ki'p  |3)  +  . . .  +  e2”v  |7))  | q) 

3.  (|0)  +  e2™*  |1)  +  e2ni2v  |2)  +  . . .  +  e2™3^  |7))  | q) 

The  possible  multiplicities  j  of  p  in  the  coefficients 
e are 


III.  QUANTUM  ALGORITHMS  WITH  POWER 
QUERIES  FOR  PHASE  ESTIMATION 

We  would  like  to  derive  lower  bounds  for  any  quantum 
algorithm  with  power  queries  that  estimates  the  phase  of 
a  matrix  Q  for  a  given  eigenvector  Ig).  In  other  words 
the  set  of  allowed  inputs  for  our  problem  is 


Ji  =  {0,  p\ ,  p2  ,  Pi  +P2}  =  (0, 1,2,3}. 

4.  (|0)  +  e2^  |1)  +  e2vi2tp  |2)  +  . . .  +  e2jli7cp  |7))  | q) 
The  possible  multiplicities  after  this  step  are 

Js  =  {j,j+P3  ■  j  €  J2}  =  {0, 1, ..  .,7}. 


Q\q),t  =  { Q  :  Q  is  a  unitary  t  qubit  transform, 

|g)  is  an  eigenvector  of  Q}. 

We  now  give  a  framework  that  is  general  enough  to 
allow  us  to  analyze  any  algorithm  with  power  queries 
WF(Q)  that  solves  this  problem.  The  most  general  algo¬ 
rithm  will  be  of  the  following  form: 

|V-(T)(0))  =  UtW£(Q)Ut-  1  • . .  W£{Q)U0\1>W).  (4) 

Here  the  Uq,  U\,  . . .,  Ut  are  arbitrary  but  hxed  c+t  qubit 
unitary  transformations  and  a  fixed  c  +  t  qubit 

state,  for  example  |0)  |g).  In  our  analysis  we  neglect  the 
cost  to  implement  the  Uj  or  to  prepare 

The  varying  parts  of  algorithm  (4)  are  the  Wf1  = 
Wi1  (Q):  power  queries  of  Q  for  pj  G  N,  lj  =  l,...,c, 
and  c  G  N  arbitrary.  A  measurement  of  the  final 
state  |V’(T)((3))  in  the  standard  basis  yields  a  state  \k), 
k  =  0, . . . ,  2c+t  —  1,  with  probability  Pk,Q ,  from  which 
we  get  a  solution  p(k)  G  [0,  1).  If  for  all  Q  G  Q\q),t  the 
probability  to  get  an  e-estimate  to  the  correct  phase  tp  of 

Q 

Pk’Q  -  i’  (5) 

k:\\ip-^(k)\\<e 

then  the  algorithm  (4)  solves  the  phase  estimation  prob¬ 
lem  to  within  e  with  probability  |  in  T  power  queries. 

We  are  interested  in  the  smallest  number  T  such  that 
a  quantum  power  query  algorithm  of  form  (4)  fulfills  con¬ 
dition  (5)  for  all  Q  G  Q\q),t- 


The  final  step,  the  inverse  Fourier  transform,  does  not 
depend  on  p.  It  also  does  not  change  the  possible  mul¬ 
tiplicities  of  p,  but  just  creates  linear  combinations  of 
them.  Consider,  e.g.,  the  coefficient  of  the  state  |2): 

r  7 

Y,e-27ri2j/8e2nijvW\q)  =  £  e2irij(<p-l/4)  |2^  , 

j=0  1=0 

which  gives  the  probability  p2{lp)  of  measuring  |2): 


P2{?)  =  |  £  g27rij(vj-l/4) 

3=0 


2 


7 

e27r*(j— Z)(v» — 1/4)  ^ 

3,1=0 


which  is  plotted  in  figure  2.  Figure  2  shows  that  the 
probability  that  |2)  is  measured  is  high  if  p  is  close  to 
0.25,  which  is  the  value  represented  by  |2). 

The  figure  indicates  that  the  width  of  this  probability 
peak  depends  on  the  frequencies  present  in  the  proba¬ 
bility  function:  higher  frequencies  allow  sharper  peaks. 
The  goal  of  this  paper  is  to  prove  that  every  halving  of 
the  width  of  the  probability  peak  requires  one  additional 
step  of  the  algorithm. 

The  proof  consists  of  three  steps.  The  first  is  to 
quantify  the  influence  of  each  additional  application  of 
Wf  on  the  frequencies  present  in  the  probability  func¬ 
tion  (Theorem  2).  Now  consider  a  probability  function 
as  in  figure  2.  It  must  have  a  high  peak  >  3/4  with 
small  width  e  and  should  be  close  to  zero  everywhere 
else.  We  will  show  that  such  a  function  requires  a  large 
range  of  frequencies  to  be  present  (lemma  3).  Finally  we 
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FIG.  2:  (Color  online)  The  probability  of  measuring  the  state 
1 2)  depending  on  ip  for  the  algorithm  depicted  in  figure  1  with 
T  =  2  and  T  =  3. 


will  show  an  upper  bound  on  the  number  of  frequencies 
that  a  power  query  algorithm  can  generate  with  T  power 
queries,  which  proves  Theorem  1. 

In  Theorem  2  consider  a  general  Q ,  which  acts  on  t 
qubits  and  therefore  has  2*  eigenvectors  | ips)  with  eigen¬ 
values  e2m,fia.  We  assume  that  the  eigenvectors  \ips)  are 
fixed,  but  that  the  eigenvalues  change.  We  will  prove 
that  after  T  steps  only  coefficients  like 

ae 2iri(ji  <Pl+—+j3t  <p2t ) 

will  occur.  Here  (j i, . . .  ,  jj)  is  from  the  set  Jr  defined 
by  the  recursion 


Jt+1  ■=  {  (jl ,  {jl  +PT+ 1,  -  -  - ,  J2*  ) ,  •••, 

+PT+l)  ■  {jl,  ■~,j2t)  G  Jt}  (6) 

and  Jo  =  {(0, . . .  ,0)}. 

Theorem  2.  Let  Q  be  a  unitary  operation  with  eigen¬ 
vectors  \ijja)  and  corresponding  eigenvalues  e'2'Kl‘Pa ,  s  = 
1 , . . . ,  2* .  Let  the  \ips)  be  fixed  and  vary  the  phases 
ips  G  [0,  1).  Any  quantum  algorithm  with  power  queries 
Wp  =  WP{Q),  fixed  unitary  transformations  Uj  and 
starting  state  can  be  written  as 

UTW%  ...UJV^U^)  =  $>£%  1.-.^')  |fc>  (7) 

k 

(T) 

for  all  ps  G  [0,  1),  where  the  Sj.  ypi ,  ■  ■  • ,  <^2‘ )  are 
trigonometric  polynomials  of  the  following  form: 


q(T) 

‘-’fe 


{Vi,- 


■■,<P2*) 


-  E 

(ji,--,32t)&JT 


AT)  2-Ki(jlVl+...+j2t<p2t) 

ak,(h,...,j2t)e 


(8) 


with  a 


(T) 

Mil 


■  6  C  and  Jt  defined  by  (6). 


To  shorten  our  proofs  we  will  use  the  following  short¬ 
hand  notations.  We  say  that  a  trigonometric  polynomial 
S^\<p i, . . . ,  (/?2‘)  is  of  powers  Jt,  indicated  by  the  ^ 
superscript,  if  it  can  be  written  as  a  sum  over  Jt  as  in  (8). 
Furthermore  we  abbreviate  the  vectors  j  =  {j i , . . . ,  jj ) , 
(p  =  (ipi, . . . ,  and  define  the  following  notation: 


j  -0  =  jl<Pl  +  ■  ■  •  + 


Proof  of  Theorem  2:  To  simplify  the  proof  we  choose  a 
different  basis  instead  of  the  standard  basis.  We  divide 
the  quantum  state  into  a  control  part  |to)  and  an  eigen¬ 
vector  part  | ips)  and  show  that  we  can  write 


UtW[tt ...W^U0\^)  =  Y,  E  Ms) ,  (9) 


m= 0  s= 1 


for  trigonometric  polynomials  of  powers  Jt, 


=  E 

ie  Jt 


g(T)  g2 niff' 

m,s,j 


(10) 


The  proof  is  by  induction  on  the  number  of  queries  T. 
For  T  =  0  we  have  no  dependence  on  (p  and  is 

just  a  constant  since  Q  was  not  applied: 


U0 


2C  — 1  2* 

=  E  E^N^)  im>^) 

m=0  s=  1 

=  :  E  E“m!s,(0 . 0)  KV’s) 

m= 0  s= 1 

m=0  s=l  JGjo 


Let  now  T  be  arbitrary  and  let  equations  (9)  and  (10) 
hold.  If  we  apply  to  (9),  only  those  states  | m,ipB) 

are  affected  for  which  the  Zt+ i-th  control  bit  is  set,  i.e. 
m;T+ 1  =  1.  For  these  states  we  get 

wW+i  lTO>  V’s)  =  |m)  QPT+1  \ips)  =  | m)  e2mpT+lVa  |^s) 
and  therefore  the  coefficient  of  | changes  to 


S^s(0)e2”ip^ 


E 

IS  Jr 


pfi,T)  ^2-Kif-ip e2-KipT+TPs 
m,s,j 


E 

j&Jt 


g(T)  ^elTti{jnpi  +  ...+(jB+PT+l)<Ps  +  ---32tV2t) 
m,  s,j 


But  this  is  a  trigonometric  polynomial  of  powers  Jr+i 
(recall  eqn.  (6)),  and  we  can  write  it  as 


s£+1]{fi) 


E 


al'T+1)  e2ni3'V 
m,s,  j 


3&Jt+ i 


4 


if  we  define  the  coefficients  c/T+-*  properly: 

m,s,j  ^  ^  J 

5(t+1)  =  a(T) 

for  pt+i  <  js  and  0  otherwise. 

Finally  we  define  <Sy^E(<^)  :=  Smtl(p)  for  the  states 
for  which  the  control  bit  is  not  set  (m/t+1  =  0)  and  we 
can  write 

2C  — 1  2t 

Kt+iUt  . . .  WKUotyW)  =  E  E  ®  1}(0  K  ^ )  • 

m— 0  s— 1 

Now  we  use  that  the  transformation  Ut+i  and  the 
eigenvectors  |'0S)  are  fixed  for  all  algorithms  we  consider 
and  define  uym  s.rlt  =  (m,  ifs\  Ut+ i  \n,ipt)-  Then 


up  to  precision  e  has  to  use  power  queries  W?* , 

WfT  such  that  the  set 

JT 

mt  =  Epfc|^c{i,...,r}}  (13) 

feeic 

has  more  than  \A4t\  >  ^  elements. 

Note  that  this  is  both  a  condition  on  pi , . . . ,  pr  as  well 
as  on  T:  by  cleverly  choosing  the  powers  pi, ...  ,pr  we 
can  get  away  with  a  smaller  T,  e.g.  for  pj  =  2J~1  we  get 

Mt=  {— 2t  +  1,  — 2t  +  2, . . . ,  2t  —  1}  , 

while  for  the  choice  pj  =  1  we  only  get 

Mt  =  {-T,-T+1,...,T-1,T}. 

Proof  of  lemma  3:  For  simplicity  assume  that 


%Es£+1)®in>w 

n,t 

=  E  Sn/1\0)UZ+sln,t  \™,1ps),  (H) 

m,s,n,t 

which  gives  the  following  coefficient  for  \m,ips)' 


2e  =  —  for  some  N  €  N.  (14) 

N 

We  will  analyze  the  behavior  of  all  possible  algorithms 
for  the  phase  estimation  problem  on  a  special  subset  of 
Q\q),f  Fix  some  arbitrary  vectors  \>p 2) ,  ■  ■  ■ ,  \ip2t)  such 
that  \q) ,  \if>2) ,  •  •  ■ ,  | Jp2t)  form  an  orthonormal  basis  and 
consider  the  following  input: 


E  E 


(T+1)  ,,(T+1) 


cr  E e 

n,t,j 


j&jT+l 


(T+i)  ~(T+1) 


02ixij-(p 


=  E  E 

JeJr+1  n't 

=■  E  a^+s%^=:S^t1](p)- 


(12) 


ie  Jt+i 


This  completes  the  induction  and  establishes  equations 
(9)  and  (10). 

Using  the  same  argumentation  as  in  equations  (11)  and 
(12)  we  can  finally  rewrite  the  state  in  equation  (9)  in  the 
standard  basis  |fc)  e  { |0) ,  |1) , . . . ,  |2c+t  —  1)}  through 


a 


(T+l) 

k,j 


E  ( k\rn ,  ips) 


~(T+ 1) 

a  4 

m,s,j 


m,s 


(rp\ 

and  Sk  (p)  is  of  the  same  powers  as  Sm,s(p),  which  is 
of  powers  Jt-  This  proves  equation  (7)  and  (8).  □ 


We  now  focus  on  the  specific  problem  of  phase  esti¬ 
mation.  The  next  lemma  provides  us  with  a  necessary 
condition  on  the  powers  pi,p2,  ■  .  ■  such  that  a  quantum 
algorithm  with  power  queries  can  solve  the  phase  estima¬ 
tion  problem  with  precision  e. 

Lemma  3.  Any  quantum  algorithm  estimating  the  phase 
p  of  an  eigenvector  \  q)  of  matrices  Q  from  the  class 


Q\q)  t  =  {Q  :  Q  is  a  unitary  t  qubit  transform , 

|gr)  is  an  eigenvector  of  Q }. 


2* 

Qr  :=  e2ni2re  \q)  <«|+E  ^iVa  E  W*  I  •  (15) 

s=2 

The  phase  p  we  are  interested  in  is  p  =  2 re  for  input  Qr. 

Since  the  difference  between  the  phases  of  the  matrices 
Qr  is  2e  and  we  require  e  correctness,  a  measurement  will 
yield  states  | k)  from  the  distinct  sets  Br: 

Br  =  |fc  :  ||2re  —  k\\  <  e  j .  (16) 

Depending  on  the  number  of  qubits  we  use  in  our  quan¬ 
tum  algorithm,  the  sets  Br  can  contain  one  or  more 
states. 

By  Theorem  2  we  know  that  we  can  write  the  coeffi¬ 
cient  of  a  state  |fc)  €  Br  after  T  queries  as 


(p,p2,...,p2t) 


^  E 


(T)  2-Ki(j1ip+j2<p2+...+j2t‘P2t) 

ak,Uu-,32t) 


In  this  proof  we  are  only  interested  in  the  behavior  for  the 
Qr.  Therefore  we  can  drop  the  dependence  on  p2,  ■  ■  ■ ,  p2 « 
and  let 

EE)  :=  S4\p,p2,...,p2t)  =  E 

IGCt 

where  Cq  =  {0}  and 

Tt  ={ji  '■  {ji,  ■  ■  ■  > 32*)  G  Jt) 

=  {j%,  ji+PT  ■  (ji,  •  •  •  ,72*)  6  Jt-i) 

-{E/'a-  A  C  { 1 , . . . ;  7  }  J 

keK 
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and  the  coefficients  (note  that  l  is  fixed  on  the  right  side) 


state  from  Br  is  correctly  returned  from  PBr(p), 


/3 


(T) 

k,l 


E 


nC1")  2t ri(j2¥>2  +  ...+j2tv2t) 


The  probability  PBr( <fi)  of  measuring  a  state  from  the 
set  Br  defined  in  (16)  of  all  e  approximations  to  ip  =  2 re 
is  now  given  by: 


PbA1?)  ■=  E  Sk  V) 


keBr 


(T)  g27ri(Z  —  l')ip 


=  E  E  E 

k(zBr  IeCt  1'€jCt 

=TE  At 
=■■  E  ^ 

mG-MT 


(^)  ft(T)  ^2Tvi(l— l')<p 


(17) 


1(T)  e 2irim<p 


with  the  set  .Mt  given  by 


AfT  =  {Z-/'  |/,f'e£r} 


and  the  coefficient 


„(T) 

Ir.m 


E 


'r,iP 


i,i’ecT 

l—l'=m 


E  EOS’- 

l,V G-Ct  k£Br 
l—l'=m 


For  illustration  recall  figure  2,  which  shows  exactly  one 
of  these  probability  functions  psr  (p)  and  also  their  highly 
oscillatory  behavior.  In  the  case  of  the  phase  estimation 
algorithm,  Br  =  {|r)|  and  the  figure  shows  Pb2  ((/?)■ 

We  apply  the  Discrete  Inverse  Fourier  Transform  to 
PBr(p)  (evaluated  at  the  points  ip  =  n/N  for  n  = 
0 , ,N  —  1)  and  get  for  the  fc-th  coefficient 


JV-l 


E«*4> 


n— 0 


e  N 


> 


PBr(^)e  « 


(  N-l 

E1  ,  Tl  .  1-Kikn 

PBr(^)e-^~ 

n—0 

n^r 


(20) 


”1  N -1 

O  v  \  ,  ll  , 

— z —  E 


n— 0 

riy£r 


since  the  probability  PBrip)  has  to  obey  PBr(r/N)  >  | 
(recall  the  definitions  of  pBr  and  Qr).  If  we  knew  that 
for  the  second  term  in  (20) 


J2pBr(n/N)  <  \ 

n— 0 
n^r 


we  could  establish  that 


(21) 


N-l 

E,  ,  r \  _ 2nikn 

Pbt  {n/N)e  N 

n—0 


N  E  ^rTm 

wiGMt 
m=k  mod  N 


>  o. 


We  will  show  that  this  property,  while  not  necessarily 
always  true,  will  be  true  for  at  least  most  of  the  pBr(p)- 
There  are  N  different  possible  outcome  sets  B0 ,  ..., 
Bn- i-  We  know  that  for  any  p  =  n/N  all  mutually 
exclusive  probabilities  of  measuring  a  state  from  B,  for 
r  =  0, . . . ,  N  —  1  have  to  add  up  to  at  most  1: 


N-l 

Y  psr{n/N )  <  1  for  n  =  0, . . . ,  N  —  1, 

r— 0 

Let  R<  C  {0, . . . ,  N  —  1}  be  the  set  of  all  r  for  which 
(21)  holds  and  R-  the  set  for  which  it  does  not.  |i?<| 
has  to  be  greater  than  1  since  we  can  split 


JV  — 1 

yPBr(n/N)e-2vikn/N 

n—0 


N-l  N-l  N-l 

N  =  E 1  ^  E  E^(^) 

n—0  n—0  r=0 


N-l 


=  E  E  Nrn6 


(T)  2ni(m—k)n/N 


n—0  mG A^t 


=n  y  ?7(t) 

1 Y  /  j  'lr,mi 


rnGJ^T 
i=k  mod  N 


(18) 


since  for  m  ^  k  mod  N 


N-l 


^  ^  ^27ri(m—k)n/N 
n= 0 


^  _  g2ni(m—k)N/N 

^  _  g27ri(m— k)l/N 


(19) 


We  can  bound  (18)  by  separating  the  part  where  a 


into  the  following  parts: 


N-l  N-l  N-l 

E^r(^)+  E  E^r(^)+  E  E^(^f) 

r— 0  rG R<  n—0  rGi?-  n—0 

n^r  n^r 

>E+o|f?<|  +  ^|^|  =  E+^|^i 

-  4  i  i  4  i  i  4  4  1  1 

and  therefore  |f?-|  <  | N  or  \R<\  >  |7V. 

Thus  there  is  an  r  €  R<  for  which  eqn.  (21)  holds  and 


0  < 

"  E 

„(t) 

'lr,m 

<n  y 

m 

'  lr,m 

mG 

wGMt 

(22) 


m—k  mod  N 


m—k  mod  N 
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for  all  k  =  0, . . . ,  N  —  1. 

(T) 

This  means  at  least  N  of  the  r/fm  have  to  be  nonzero 
and  thus  PBr(<p)  from  eqn.  (17)  must  have  at  least  N 
nonzero  terms.  In  other  words 

|A4T|>iV=^,  (23) 

where  we  used  the  definition  of  N  in  (14).  □ 

Lemma  3  now  allows  us  to  give  a  lower  bound  for  the 
phase  estimation  problem.  The  numbers  of  queries  T 
that  is  needed  by  any  quantum  algorithm  is  D(log  -  ). 

Proof  of  Theorem  1:  From  lemma  3  we  know  that  the  set 
Mr  has  to  have  \Mt\  >  ^  elements.  We  can  easily 
derive  an  upper  bound  on  |A4t|-  There  are  at  most  2T 
elements  in  the  set 

C-t  =  {  Y,Pk\KC{l,...,T}} 
keK 

and  therefore  we  have  |A4t|  <  |£t|1 2 3 4 5 6 7 8 9  <  (2T)2  =  22T. 
Combining  our  estimates  for  |AIt|  we  get 

22T  >  \Mt\  > 

and  the  number  of  queries  T  must  grow  like 
T>  ^log2^  =  fi(log^). 

□ 


V.  CONCLUSIONS  AND  EXTENSIONS 

In  this  paper  we  have  obtained  lower  bounds  for  quan¬ 
tum  algorithms  that  approximate  the  phase  estimation 
problem  through  the  new  lower  bound  proof  technique 
of  frequency  analysis.  These  lower  bounds  match  the 
known  upper  bounds  for  phase  estimation. 

The  frequency  analysis  method  can  be  used  to  give 
lower  bounds  for  the  Sturm-Liouville  eigenvalue  problem 
[7].  The  results  also  extend  to  other  forms  of  quantum 
queries,  e.g.  the  query 

Qf  |z>  I V)  =  \x)  |  V  ©  f{x))  =  \x)  | y  +  f{x)  mod  2m) . 

An  application  of  the  frequency  analysis  method  to  these 
problems  will  be  the  subject  of  future  work. 
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